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Abstract. Accurate molecular spectroscopy in the mid-infrared region allows preci¬ 
sion measurements of fundamental constants. For instance, measuring the linewidth of 
an isolated Doppler-broadened absorption line of ammonia around 10 jim enables a de¬ 
termination of the Boltzmann constant ^b- We report on our latest measurements. By fit¬ 
ting this lineshape to several models which include Dieke narrowing or speed-dependent 
collisional effects, we find that a determination of with an uncertainty of a few ppm 
is reachable. This is comparable to the best current uncertainty obtained using acoustic 
methods and would make a significant contribution to any new value of determined 
by the CODATA. Furthermore, having multiple independent measurements at these ac¬ 
curacies opens the possibility of defining the kelvin by fixing an exciting prospect 
considering the upcoming redefinition of the International System of Units. 


1 Introduction 


Atomic physics has brought us some of the most precise measurements of fundamental constants 
ina, stringent constraints on their variations in time lO, clocks of the highest accuracy lH, quantum 
sensors II5I6L probes of general relativity 13 and physics beyond the standard model of particles and 
fields 0. This formidable discipline has historically gone hand in hand with the study of molecules, 
which contributed to the development of laser physics and high-resolution laser spectroscopy, showing 
for instance the first observation of the photon recoil splitting 0. Molecular systems provided the first 
natural-resonance-based frequency standard in 1949 using microwave transitions in ammonia ifTOl . 
Optical frequency chains developed to measure molecular transition frequencies led to the accurate 
measurement of a fundamental constant, the speed of light ifTTIl . Note in passing that the need for higher 
resolution and frequency accuracy inspired the first ideas on the use of cold atoms and molecules 
lfT2TT3 . The potential for molecules to outperform precision measurements carried out on atoms, 
owing to their rich electronic, vibrational, rotational and hyperfine structure is becoming increasingly 
apparent. Molecules are now being (or have recently been) used to test fundamental symmetries such 
as parity II14I15I16II and time reversal (in ifTTll . the upper limit on the size of the electron electric 
dipole moment - a signature of parity and time reversal if non-zero - was set using YbF molecules, 
outperforming the previous limit set by measurements on atoms), to test the symmetrization postulate 
of quantum mechanics ifTSll . to measure either absolute values of fundamental constants (electron-to- 
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proton mass ratio in or Boltzmann constant as discussed in this paper), or to measure their variation 
in time (fine structure constant a ll20ll . electron-to-proton mass ratio II21I22II ). 

The current International System of Units (SI) is built on seven well-defined base units: the meter, 
the kilogram, the second, the kelvin, the ampere, the mole and the candela. The definition of some 
of these units is based on macroscopic objects or macroscopic properties of matter. Since 1889 the 
unit of mass has been defined by the international prototype of the kilogram, a platinum cylinder 
located at Bureau International des Poids et Mesures (BIPM) in Sevres (France). This definition has 
a major drawback. It is based on the mass of an artifact, whose composition inevitably fiuctuates in 
time. In more than a hundred years, comparisons of 7C with other prototypes showed drifts on the 
order of 50 yug (5 x 10“^ in relative value). This directly impacts the definition of the mole, the ampere 
and the candela which depend on it. Just as unsatisfactory, the kelvin is defined by the temperature 
(273.16 K) of the triple point of water of a specific isotopic composition. The precision with which 
this isotopic composition can be achieved currently limits the reproducibility of the kelvin realization 
at a few 10“^. A revision of the SI is thus expected in the near future II23I24I25I26I27II . One objective 
is to avoid using artifacts or macroscopic properties of matter. Another objective is to redefine units 
in terms of fundamental constants expected to not vary in time. A third objective is to derive base 
units from the time unit which is by far realized with the best accuracy. This was done in 1983 for 
the length unit. The velocity of light c was fixed to 299 792 458 m • s“\ defining thereby the metre, 
as the length travelled by light in vacuum in 1/299 792 458 s. In the new SI four of the base units, 
namely the kilogram, the ampere, the mole and the kelvin, will be redefined, based on fixed numerical 
values of the Planck constant h, the elementary charge e, the Avogadro constant Na and the Boltzmann 
constant ks, respectively. Let us detail the kelvin case. Temperature has a microscopic interpretation. 
It is proportional to the mean energy per particle per degree of freedom E = ^k^T. hjE having the 
dimensions of time, fixing the value of h and is a way to define temperature by linking it to the time 
unit. The revision of the SI has been envisaged thanks to remarkable progress made over the last few 
years towards the experimental determination of /z, Aa and ^b- 

We report on our latest progress on the determination of the Boltzmann constant ^b by accurate 
molecular spectroscopy in the mid-infrared region, using the Doppler Broadening Technique (DBT). 
The 2010 value recommended by the CODATA (Committee on Data for Sciences and Technology) is 
^B( 20 i 0 ) = 1.380 648 8(13) x 10“^^ J.K“\ with an uncertainty of 0.9 ppm, based on eight published 
values ll^ . Six of them were obtained using a single method, namely the Acoustic Gas Thermometry. 
We are aiming to provide a competitive measurement using an alternative method in order to contribute 
to the next CODATA value and to the new definition of the kelvin. 


2 The Doppler Broadening Technique 

2.1 Principle 

The principle of the DBT 0, illustrated in Figure [T] is to record the Doppler profile of an absorption 
line of a vapour in thermal equilibrium. In the Doppler regime, when all other broadening mecha¬ 
nisms are negligible, the recorded Gaussian profile derives from the Maxwell-Boltzmann distribution 
of velocities along the laser beam axis. The Boltzmamn constant is then calculated fromzlvo the half- 
Doppler width at 1 / ^ : 

/^^d\ 
h ~ hT\v^] 

where m is the molecular mass, vq the central frequency of the molecular line and T the gas tempera¬ 
ture. 

Using saturated absorption spectroscopy, we measure vq with a relative uncertainty of a few 10“^ 
ll29l . The Watt balance experiment measures h at the level of 5 x 10“^ ll^ and him has an uncertainty 
of 10“^ deduced from atom interferometry experiments with rubidium O and atomic mass ratios 
measured in ion traps m . The uncertainty on ^b is then dominated by uncertainties on the temperature 
and the Doppler width. The key ingredients required for an accurate determination of are the 
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thermostat at 273.15 K 



following: (i) control and accurate measurement in real time of the gas temperature, by placing the 
absorption cell in a thermostat, and (ii) a precise record and reliable model of the lineshape in order 
to extract zIvd from the data. In Section [3] we discuss how we are able to control and measure the 
temperature. Section |4| details the Doppler width measurement and its accuracy. 


2.2 The experimental setup 

The experimental setup has been described in detail in II31I32II . It contains a CO 2 laser stabilized on 
a saturated absorption line of osmium tetroxyde. Thanks to the resulting unique spectral properties 
(below 10 Hz width, 0.1 Hz frequency instability for 100 s integration time O), an incomparable 
control of the frequency scale is reached, a crucial element for a careful reproduction of the lineshape. 
The intensity of the laser beam is stabilized before being sent to the absorption cell filled with gaseous 
ammonia. NH 3 was chosen because of its strong V 2 absorption band which overlaps with the 8-12 yum 
spectral range of our ultra-stable spectrometer. Furthermore, this band exhibits well-isolated transi¬ 
tions which minimize any deformation of the line shape due to neighbouring lines. The absorption 
length of the cell can be adjusted from 37 cm in a single-pass configuration (SPC) to 3.5 m in a mul¬ 
tipass configuration (MFC). By varying the pressure from 1 to 25 Pa (SPC) and 0.1 to 2.5 Pa (MFC), 
absorption amplitude can be varied between 20 to 98%. The laser frequency is tuned close to the y 2 
saQ(6,3ii| line of ^“^NHs (of measured central frequency vq = 28 953 693.9(1) MHz), and scanned 
via a tunable electro-optic modulator. The absorption cell is placed inside a thermostat at 273.15 K 
(see Section O. In order to avoid pressure induced modification of the lineshape (see Section 14.1.11) 
while preserving a good signal-to-noise ratio, spectra are recorded at pressures around 1 Pa in MFC 
(the SPC is used when pressure effects are to be studied). Typical spectra are displayed on Figure [T] 
They are recorded in steps of 500 kHz with 30 ms of integration time per point, leading to a typical 
signal-to-noise ratio of 10 ^. 


3 Temperature control 

This experiment requires the molecular gas to be maintained at a constant and homogeneous temper¬ 
ature. A thermostat has been designed in close collaboration with thermometry experts from Labo- 

^ J = 6 and K = are, respectively, the quantum numbers associated with the total orbital angular momentum 
and its projection on the molecular symmetry axis. 
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ratoire Commun de Metrologie, LNE-CNAM. The absorption cell is placed inside a copper thermal 
shield which is itself inside a stainless steel enclosure immersed in about 1 m^ of an ice-water mixture 
which stabilizes the temperature close to 273.15 K. The thermal shield is only linked to its enclosure 
via a copper ring. The absorption cell is mechanically linked to the thermal shield via Teflon hold¬ 
ers to ensure that heat transfer between them is predominantly radiative. This, in turn ensures that 
the heat transfer is spatially smoothed, thereby reducing residual temperature gradients, but also in¬ 
creases the thermal time constant of the system (to about 10 h), which smoothes temperature drifts 
over similar time scales. The gas temperature is measured with Hart capsule standard platinum resis¬ 
tance thermometers (cSPRTs) directly placed inside vacuum, in contact with the absorption cell. Those 
cSPRTs are calibrated at the triple point of water and at the gallium melting point, and are coupled to 
a resistance measuring bridge (Guildline Instruments Limited, Ref. 6675/A) calibrated against a resis¬ 
tance standard with a very low-temperature dependence. In this thermostat the temperature of the gas 
probed by the laser beam shows an overall accuracy of ~ 1 ppm, deduced from a thorough analysis of 
the thermal stability and homogeneity of the global temperature control system ll^ . 


4 Doppler width measurement 

The Doppler width is obtained by fitting recorded spectra (see Eigure[T]) to the Beer-Lambert law (with 
an added baseline): 

P(y) = [^0 + s(v - Vo)] exp[-A/(y - vq)] (1) 

with Pq the incident power, 5* the slope of the baseline, v the laser frequency, A the integrated ab¬ 
sorbance (proportional to the pressure and the absorption length) and I(y - vq) the normalized ab¬ 
sorption profile that contains the Doppler width Avd. The exact form of the profile depends on the 
assumption made for the type of collisions between molecules, as will be discussed in Section IT. 1.11 


4.1 Accuracy 

We have listed and investigated a number of systematic effects that may affect our determination of the 
Boltzmann constant. Systematic shifts in our measurement have to be measured or estimated in order 
to either cancel or correct for them. To build a proper error budget, it is also essential to estimate the 
uncertainty on each correction. 


4.1.1 The lineshape problem 

When the pressure goes to 0 Pa, the collision rate goes to 0. The absorption profile I(y - vq) (see 
Equation[T]) is only Doppler broadened and is a simple Gaussian profile of half-width at 1/e Avd, of 
about 50 MHz in our experimental condition^ In practice however, collisions induce both an addi¬ 
tional broadening and a shift. The overall lineshape is then a Voigt profile (VP), i.e. the convolution 
of the above mentioned Gaussian profile with a Lorentzian profile. In our experimental conditions the 
collisional broadening and shift (the actual half-width at half-maximum y and shift S of the Lorentzian 
function) are typically 100 kHz/Pa and 1 kHz/Pa respectively. The VP is probably the most widely 
used model in the spectroscopy community. However more subtle effects cannot be ignored in our 
case. Collisions induce an auto-confinement of molecules by the surrounding molecules which leads 
to a narrowing of the Doppler contribution to the lineshape . This narrowing is important when 
the pressure is such that this confinement is at the scale of the wavelength, the so-called Lamb-Dicke- 
Mossbauer (EDM) regime. In our case the EDM narrowing is expected to be about yS ~ 10 kHz/Pa 
El, where p is the frequency of velocity-changing collisions, an additional parameter of the line- 
shape. Depending on the assumption made for molecular collisions, adding this contribution to a VP 
leads to different lineshapes. The Galatry profile (GP) assumes soft collisions against light particles 

^ typical natural width for rovibrational levels are negligible, of the order of a few Hz. 
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Fig. 2. Normalized residuals for nonlinear least-squares fits of the same average of 24 spectra, recorded at 17.3 Pa, 
to the exponential of a Voigt, a speed-dependent Voigt, a Galatry and a Rautian profile. 


in a typical Brownian motion picture, in which the memory of the initial velocity is conserved after 
many collisions 1^ . The Rautian profile (RP) assumes strong hard sphere collisions against infinitely 
massive particles, after each of which the memory of the initial velocity is lost II37I38L Another subtle 
effect is the fact that the collisional broadening and shifting actually slightly depend on the molecular 
speed. Taking this into account leads to the speed-dependent Voigt profile (SDVP) with an overall 
narrowing compared to the usual Voigt profile. The impact of all collisional effects considered here 
are small in the typical low pressure experimental conditions favored for the determination of . In 
the general case, the lineshape can result from the combination of all these effects leading to more 
complex lineshapes 1 ^ . 

Knowing the true lineshape of a given absorption line is fundamentally difficult. Until recently, 
the main source of uncertainty on determining was the lack of knowledge of the correct collisional 
model for the gas, and therefore the correct lineshape to use for fitting to our experimental data. It 
for instance, constituted 97% of the 144 ppm standard uncertainty published in ll29ll . The influence 
of the lineshape model on the spectroscopic determination of the Boltzmann constant has been raised 
by us and other groups working with O 2 and H 2 O II29I32I40I4 1 142I43I44I45II . We recently carried out 
an extensive lineshape study at relatively high pressure in SPC, in order to enhance collisional effects 
and to determine the most suitable model for our data ll40ll . The residuals obtained after fitting the 
average of 24 spectra recorded at 17.3 Pa to the exponential of a VP, SDVP, GP and RP are displayed 
in Figure [2j Residuals for the VP, GP, and RP exhibit structures with amplitudes respectively 10.3, 
4.3 and 6.7 times larger than the rms noise level, clearly indicating that these profiles do not match 
the experimental data. By contrast, analysis with the SDVP demonstrates a good agreement between 
this model and the experimental lineshape. In the covered pressure range, the deviation from a Voigt 
profile of the y 2 saQ(6,3) line of seems dominated by speed-dependent effects and not LDM 

effects. The following speed-dependent lineshape parameters have been obtained at the few percent 
level: y = 120(3) kHz/Pa, 6 = 1.2(1) kHz/Pa, m = 0.360(9) and n = -3.8(30 


^ the m and n parameters are related to the speed-dependence of the collisional broadening and shift: under the 
hypothesis that the molecular interaction potential V is proportional to some inverse power of the internuclear 
separation R, V ~ l/R*^, where q is determined by the kind of the interaction, the dependence of the collisional 
width and shift versus the relative collisional speed are given by r{v^) ~ and A{v^) ~ respectively ESI. 
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We have clearly identified the SDVP as the most suitable lineshape model. Having extracted the 
lineshape parameters m and n with good accuracy from data taken at high pressure, we fix these param¬ 
eters when fitting the low pressure data used to determine • The uncertainty on due to collisional 
effects is then only limited by our knowledge on the associated parameters m and n. Simulated SDVP 
spectra corresponding to typical experimental conditions for ^b determination have been fitted taking 
into account uncertainties on m and n (see above) in order to estimate the impact of these uncertain¬ 
ties on the Doppler width accuracy. We followed the line-absorbance based analysis method recently 
proposed in (461, which leads to an estimated uncertainty of 0.9 ppm (respectively 1.8 ppm) on zIvd 
( respectively ^b), if spectra are recorded at ammonia pressures below 2 Pa. 


4.1.2 Hyperfine structure 

The saQ(6,3) line has been chosen because it is a well-isolated rovibrational line with long-lived 
levels. However, owing to the non-zero spin values of the N and H nuclei, an unresolved hyperfine 
structure (HFS) is present in the Doppler profile. It causes a broadening of the lineshape according to 
the relative position and strength of the 78 individual hyperfine components. Although the 12 strongest 
components are spread about 50 kHz, weaker lines around + ~ 600 kHz away from the main structure 
must also be considered, as they give the largest contribution to the HFS induced broadening. The 
relative increase in the linewidth scales as the square of the ratio zlvhyp/zIvD, with zlvhyp being the 
spread of the HFS. For zIvd ~ 50 MHz, the impact is a few ppm. 

Accurate knowledge of the HFS is required in order to apply the right correction to the Doppler 
width. Ultra-high resolution saturated absorption spectra of the saQ(6,3) line were measured in a 3-m 
long Fabry-Perot cavity (291 . A combined analysis of our infrared spectra and microwave data from 
the literature led to an accurate determination of the hyperfine constants of the transition. The position 
of the 78 components and their relative intensities were then calculated, together with the associated 
uncertainties. In order to estimate this HFS-induced systematic effect, we fit to a single SDVP a sim¬ 
ulated spectrum composed of the sum of 78 SDVP with positions and intensities corresponding to the 
precisely determined HFS. This leads to an overestimation of zIvd of 4.356(13) ppm. The correction 
to be applied on the value of ^b is thus -8.71(3) ppm, with negligible uncertainty. 


4.1.3 Other systematic effects 

A number of other systematic effects have been considered and shown to affect the error budget of our 
measurement at a negligible level (29l . Laser-power-related systematic effects such as non-linearity 
in the detection setup, saturation of the probed transition and differential saturation of the different 
hyperfine components were simulated and investigated experimentally. A systematic shift results from 
amplitude modulation of the laser beam at 40 kHz used for noise filtering. It can be readily estimated 
from simulations with negligible uncertainty. Typical lifetimes for rovibrational levels are of the order 
of 100 ms, translating into a minute natural width. Our fitting procedure (40l is such that presence of 
outgassing or partial pressure of impurities has rigourously no impact. The world record frequency 
stability and accuracy of our CO 2 laser (see Sectionlead to negligible contributions from the laser 
linewidth or the linearity and accuracy of the laser frequency scale. 

It is sometimes mentioned that there should be an additional broadening of the lineshape due 
to the finite transit time of particles through the laser beam. It has recently been shown that this is 
strikingly not the case in linear absorption spectroscopy, provided that the medium is uniform and 
isotropic Ea. Furthermore, this absence of transit time broadening holds true even after introducing a 
speed-dependent broadening and shift in the calculations (48l . All transit effects are already included 
in the inhomogeneous Doppler broadening, independently of the optical quality of the laser beam 
wavefronts and effective diameters. Rigourously no systematic effect due to the laser beam geometry 
is thus expected on the Doppler width measurement. 
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Fig. 3. Improvement in the statistical uncertainty of measured at Laboratoire de Physique des Lasers from 2004 
onwards. 


Table 1. International state-of-the-art for the statistical uncertainty (except for Hu et al for which only the global 
standard uncertainty can be found in the literature) obtained on using the DBT. 


place 

group 

species 

wavelength {ji m) 

fractional uncertainty (ppm) 

USTC, China 

S.M. BAxetal |55l 

‘"C 2 H 2 

0.8 

2000 

NMIJ, Japan 

T. Shimizu et al ||53l 


1.5 

1200 

UWA, Australia 

A. Luiten et al ll54l 


0.8 

398 

UniNa2, Italy 

L. Gianfrani et al 15^ 

'2C'«02 

2 

90 

UniNa2, Italy 

L. Gianfrani et al 14^ 

0 

00 

1.4 

80 

UP13, France 

our group 


10.3 

6.4 


4.2 Statistical uncertainty 

7171 spectra, corresponding to 70 h of integration time, were recorded for pressures ranging between 
0.2 and 2.5 Pa. However, these were recorded in 2010, prior to the extensive lineshape study described 
in Section 14.1.11 under non-optimal experimental conditions. Regardless, the data was fitted with the 
exponential of a Galatry profile assuming the observed deviation from a Voigt profile to be a con¬ 
sequence of the LDM narrowing 1^ . The purpose of this fitting exercise was two-fold. Firstly, the 
discrepancy between zIvd values obtained by fitting Galatry and Voigt profiles estimates the systematic 
error which would occur from not knowing the true lineshape. We indeed find large systematic shift 
and a standard uncertainty of 144 ppm 1291 . highlighting the importance of seeking the correct profile. 
Secondly, it demonstrates our ability to reach a statistical uncertainty below 10 ppm. The mean of the 
7171 extracted Doppler widths was calculated to be zIvd = 49.885 90(16) MHz, leading to a statistical 
uncertainty on of 6.4 ppm. 

The improvement in the statistical uncertainty of obtained at Laboratoire de Physique des Lasers 
from 2004 onwards II29I31I32I49I50I51II is illustrated in Figure [3] To our knowledge, 6.4 ppm repre¬ 
sents the best statistical uncertainty ever obtained for the measurement of ^b using the DBT. 

Other groups followed our initiative in measuring ^b using the DBT, but probing other species 
such as CO 2 , H 2 O, C 2 H 2 and Rb II43I52I53I54I55L Table[T]gives the corresponding state-of-art for the 
statistical uncertainty obtained on ^b • Our group is currently holding the record and leading the field. 
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Fig. 4. Comparison between the eight measurements accepted to contribute to the current CODATA value of 
ED (•: Acoustic Gas Thermometry, ■; Refractive Index Gas Thermometry, ♦: Johnson Noise Thermometry). 
Vertical axis: relative difference between measurement and the current CODATA value of k^. Our next measure¬ 
ment using the Doppler Broadening Technique will have a sub-10 ppm competitive uncertainty. 


5 Conclusion 

Taking into account all systematic effects, including temperature control, leads to an overall correc¬ 
tion to be applied to with an uncertainty on this correction at the level of a few ppm. This is valid 
provided that data is recorded under the optimized experimental conditions determined by these stud¬ 
ies and provided that spectra are fitted to the SDVP, identified as the most suitable lineshape for our 
measurements. We are able to reach a statistical uncertainty on k^ of 6.4 ppm after 70 h of acquisi¬ 
tion, as demonstrated during a previous data run. Our results are considered by the 2010 CODATA 
Recommended Values of the Fundamental Physical Constants ll^ to be ”of interest” because they are 
’’obtained from a relatively new method that could yield a value with a competitive uncertainty in the 
future”. To conclude, the first determination of k^ using the DBT with a combined uncertainty below 
10 ppm is well within reach. A new series of data will soon be recorded. 

The eight measurements accepted to contribute to the current CODATA value of k^ are displayed 
in Figureldl ll^ . None of them use the DBT. Moreover, six of them were obtained using Acoustic Gas 
Thermometry (AGT). Refractive Index Gas Thermometry (RIGT) and Johnson Noise Thermometry 
(JNT) led to the two other measurements. JNT is based on the temperature dependence of the mean 
square noise voltage, developed in a resistor by thermal motion of the charge carriers. The principle of 
the gas-based measurements, AGT and RIGT, is to determine the gas-particle’s microscopic energy, 
k^T, by measuring a macroscopic temperature-dependent gas property at known pressures: the speed 
of sound for AGT or the refractive index for RIGT. In comparison, the DBT has the advantage of being 
a direct measurement of the microscopic thermal energy. Moreover, the probed atoms or molecules 
belong to a single quantum level of a well-defined isotopic species, which avoids uncertainties coming 
from macroscopic quantities, an issue for other techniques. We would also like to stress the broad and 
generic nature of the DBT, as it can be used with any molecular as well as atomic (such as ^^Rb in 
Australia) species (see Table [B. It is worth mentioning that six of the eight above mentioned values 
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were obtained between 2006 and 2010, illustrating the efforts made by several metrology institutes in 
preparation for the upcoming redefinition of the kelvin. 

A measurement using the DBT at the level of or below 10 ppm will provide another value for 
with competitive uncertainty, using an independent method. It should contribute to the next CODATA 
value and to the new definition of the kelvin. 
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